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ABSTRACT 


The  feasibility  and  accuracy  of  the  Moird  method  for 
determining  moments  and  deflections  for  laterally  loaded 
thin  flat  plates  using  Perspex  models  have  been  well  esta¬ 
blished  by  Ligtenberg  and  Bradley.  Metal  models  are  known 
to  have  some  advantages.  In  this  study,  the  Moird  method 
was  applied  to  metal  models,  e.g.  brass  and  aluminum. 

Three  cases  of  bending  problems  were  studied:  1.  sim¬ 
ply  supported  square  plate  with  concentrated  centre  load, 2. 
clamped  square  plate  with  concentrated  centre  load,  and  3. 
cantilever  plate  with  concentrated  load  at  the  midpoint  of 
the  free  edge  opposite  to  the  clamped  edge.  Theoretical  so¬ 
lutions  are  available  for  the  first  two  cases  and  the  ex¬ 
perimental  results  by  the  Moird  method  were  compared  with 
those  obtained  by  the  well  know?n  exact  solutions.  But  an 
exact  solution  is  not  available  for  the  last  case,  and  the 
results  from  finite  difference  solution  were  taken  for  the 
comparison . 

It  is  concluded  that  the  metal  models  are  applicable 
for  the  Moird  method.  This  method  provides  the  approximate 
solutions  for  the  laterally  loaded  thin  flat  plate  problems 
which  are  intractable  by  analytical  methods. 
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NOTATIONS 


Side  dimension  of  square  plate. 

Distance  from  plate  model  to  camera  lens. 

Distance  from  horizontal  axis  of  Moire  set-up  to  a  given 
point  on  the  model. 

Spacing  of  lines  on  the  grid. 

=  Et  /12(l-\±  )  ,  Flexural  rigidity  of  a  plate. 

Flexural  rigidity  for  the  calibration  plate. 

Modulus  of  elasticity 

Horizontal  distance  measured  along  the  camera  axis,  from 
plate  model  to  the  Moire  screen. 

M  Bending  moments  per  unit  length  of  sections  of  a  plate 
^  perpendicular  to  x  and  y  axes,  respectively. 

Twisting  moment  per  unit  length  of  section  of  a  plate 
perpendicular  to  x  axis. 

Concentrated  load. 

Q  Shearing  forces  parallel  to  z  axis  per  unit  length  of 
^  sections  of  a  plate  perpendicular  to  x  and  y  axes, 
respectively . 

Radius  of  calibration  plate. 

Concentrated  corner  force. 

Thickness  of  the  plate 
Deflection  of  the  plate 

z  Rectangular  coordinates. 

* 

Poisson's  ratio. 

Spacing  of  grid  used  in  the  difference  procedure. 

Slope  of  plate  model,  or  change  in  slope  at  a  point,  in 
radians . 
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I  .  INTRODUCTION 

The  general  theory  of  bending  of  a  laterally  loaded  thin 
plate  has  been  developed  for  a  number  of  special  cases.  In 
many  practical  design  problems,  owing  to  the  difficulty  of 
formulation  of  the  boundary  conditions  and  the  solutions  of 
the  partial  differential  equation  for  plates,  analytical  solu¬ 
tions  are  difficult  to  obtain.  Various  numerical  methods  or 
experimental  investigations  of  models  seem  to  involve  too  much 
work  for  solving  practical  problems  in  a  reasonably  short  time 
In  recent  years,  an  experimental  method,  the  so  called  Moire 
method,  has  been  developed.  This  method,  described  by  Ligten- 

•k 

berg  (1),  provides  sufficient  data  for  the  determination  of 
moments  in  a  plate  of  arbitrary  form  and  loading.  The  same 
method  was  also  applied  to  square  clamped  plates  with  or  with¬ 
out  square  cutouts  for  the  determination  of  moments  and  deflec 
tions  by  Bradley  (2)  (3) ,  to  flat  slabs  for  the  determination 

of  moments  by  Vreedenbugh  (15) ,  and  to  cantilever  plates  to 
determine  stresses  by  Palmer (4) .  The  application  of  Moire 
method  to  the  determination  of  displacement  was  described  by 
Weller  and  Shepard  (11) .  The  same  method  was  applied  to  the 
strain  analysis  by  Morse,  Durelli,  and  Sciammarella  (12),  and 
similar  work  was  also  introduced  by  Low  and  Bray  (13) .  The 
application  of  the  Moire  method  to  determine  transient  stress 

•k 

Numbers  in  parentheses  refer  to  the  bibliography  at  the  end 
of  this  thesis . 
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and  strain  problem  was  recently  discussed  by  Riley  and  Durelli 
(14) .  In  this  study,  the  Moire  method  was  used  to  determine 
the  bending  moments  and  twisting  moments  in  a  simply  supported 
square  plate  and  a  clamped  square  plate  subjected  to  central 
concentrated  loads.  Bending  moments  and  deflections  in  a  can¬ 
tilever  plate  subjected  to  a  concentrated  load  at  the  middle 
of  the  free  edge  opposite  to  the  clamped  edge  was  also  studied 
by  this  method. 

Results  from  this  experimental  work  were  compared  with 
Navier  solution  (6)  for  the  simply  supported  square  plate, 
Timoshenko's  solution  (6)  for  the  clamped  plate,  and  Holl's 
finite  difference  solution  (10)  for  the  cantilever  plate  res¬ 
pect  ively . 

According  to  Bradley's  work  (3)  on  the  bending  of  plastic 
plates  with  uniform  or  concentrated  loading  for  various  bound¬ 
ary  conditions,  the  agreement  of  deflection  and  moments  as 
found  by  the  finite  difference  method  and  Moire  method  was  good, 
the  maximum  variation  was  about  ten  per  cent.  Ligtenberg  found 
out  in  most  cases  (except  very  near  the  built-in  edges  of  the 
model) ,  the  accuracy  was  about  5  per  cent  of  the  maximum  moment. 
Observing  these  results,  it  is  believed  that  the  Moire  method 
is  quite  a  useful  tool  for  solving  plate  problems. 

The  aims  of  this  study  can  be  summarized  as  follows: 

a.  To  study  the  applicability  of  the  Moire  method  to  metal 
plates,  particularly  those  with  simple  or  clamped  sup¬ 
ports,  in  order  to  determine  its  accuracy. 
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b.  To  determine  the  distribution  of  moments  and  compare 
the  results  with  those  obtained  by  theoretical  methods. 

c.  To  develop  the  technique  of  this  method  by  using  metal 
models . 

t 

Although  the  basic  principles  of  the  Moire  method  were 
given  in  Ligtenberg's  article  as  well  as  Bradley's,  they  are 
reviewed  here  for  the  completeness  of  the  presentation. 

The  equipment  used  and  the  methods  followed  in  this  inves¬ 
tigation  will  be  described  in  some  detail  in  the  following 
chapters . 
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II.  THEORETICAL  ANALYSIS 

1.  Plate  Theory. 

The  ordinary  theory  of  bending  of  thin  elastic,  homogeneous 
isotropic  plates  requires  that  the  deflected  neutral  surface  of 
the  plate  shall  satisfy  the  plate  equation  (6) , 
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d  w 


-.4 

.a  w 


c)x4  +  ^c)x^c)y2  +  cVy4 
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in  which  w  is  the  deflection  of  the  plate,  q  is  the  intensity 
of  load,  and  D  is  the  flexural  rigidity  of  plate. 

The  derivation  of  equation  1  is  based  upon  the  following 
assumptions : 

a.  The  plate  material  is  elastic,  homogeneous,  and  iso¬ 
tropic  . 

b.  The  material  obeys  Hooke's  Law. 

c.  The  normal  stresses  in  the  direction  perpendicular  to 
the  plane  of  the  plate  can  be  disregarded. 

d.  Lines  perpendicular  to  neutral  surface  before  bending 
remains  so  after  bending. 

e.  Plate  must  be  in  equilibrium  at  all  times. 

f.  The  deflections  of  a  plate  are  small  in  comparison  with 
its  thickness. 

g.  There  are  no  membrane  forces  in  the  middle  plane  of  the 
plate.  This  plane  remains  neutral  during  bending. 

Now  consider  an  element  of  the  plate  which  is  bent  under 
transverse  loading.  Taking  the  positive  directions  as  shown 
in  Fig.  1,  in  which,  the  middle  plane  is  taken  to  coincide  with 
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xy  plane,  it  can  be  shown  that  (6): 
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In  these  expressions  M^.,  are  the  bending  moments  pe 
unit  length  of  sections  of  a  plate  perpendicular  to  x  and  y 
axes  respectively,  Mx^,  and  m  being  the  twisting  moments; 


6  . 


0  ,  Q  are  the  shear  forces;  V  and  V  are  the  Kirchhoff  shears, 
x  y  x  y 

and  R  is  defined  as  the  concentrated  corner  force.  Thus,  a 
complete  analysis  of  a  plate  bent  by  a  lateral  load  q  is  re¬ 
duced  to  the  integration  of  Eq.  1.  The  degree  of  difficulty 
in  performing  this  integration  depends  mostly  on  the  boundary 
conditions  existing  along  the  edges  of  the  plate. 

2.  Boundary  Conditions. 

The  edge  of  a  plate  may  be  free,  simply  supported,  elasti¬ 
cally  supported,  clamped,  etc.  In  this  study,  the  simply  support 
ed  and  clamped  edges  were  investigated.  In  the  case  of  the  canti 
lever  plate,  a  condition  of  free  edge  was  involved.  The  study 
of  these  three  boundary  conditions  will  now  be  considered, 
a.  Simply  supported  edge.  If  the  edge  of  the  plate 
along  the  y  axis,  at  x  =  0 ,  is  simply  supported, 
this  edge  is  free  to  rotate  about  the  y  axis,  and 
the  bending  moment  will  be  zero.  All  the  deflect¬ 
ions  on  this  edge  will  also  be  zero.  Therefore  at  x  =  0 


b . 


or 


x=0 


x=0 
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and  w  ~  =  0 

L  -»x=0 

Clamped  edge.  If  the  edge  of  a  plate  is  clamped, 

the  deflection  and  slope  along  the  edge  are  zero.  Let 

the  y  axis  be  clamped,  then  the  boundary  conditions  are 

U  n  =  o  \M 

L  Jx=0  loxJx_0 


=  0 


c.  Free  edge.  If  an  edge  of  a  plate,  say,  the  edge  x  =  o 
is  entirely  free.  The  boundary  conditions  existing  at 
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a  free  edge  are: 


and 
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In  order  that  a  check  of  the  experimental  results  may  be 
given,  Navier  solution  for  the  simply  supported  square  plate, 
Young's  solution  for  the  clamped  square  plate,  and  Holl's  solution 
for  the  cantilever  plate  by  the  finite  difference  method  were 
considered. 

3.  Solutions. 

Simply  supported  square  plate.  The  deflection  of  a  simply 
supported  square  plate  with  central  concentrated  load  was  given 
by  Navier  in  the  form  of  a  double  trigonometric  series  (6) , 
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*Co-ordinates  to  be  used  in  computation. 
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which  satisfies  the  plate  equation  and  boundary  conditions. 


The  series  converges  rapidly,  and  the  deflection  at  any  point 


of  the  plate  can  be  obtained  with  sufficient  accuracy  by  taking 


only  the  first  few  terms  of  the  series.  For  example,  if  the  first 


four  terms  of  the  series  were  taken,  the  maximum  deflection  is 


about  3  1/2  percent  less  than  the  correct  value. 


The  expressions  for  bending  and  twisting  moments  can  be 


obtained  from  the  general  solution  1C  by  using  Equations  2, 


3,  and  4.  These  are  as  follows: 
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where  P  is  the  concentrated  load,  \±  is  the  Poisson's  ratio,  a  is  the 


length  of  the  plate. 


As  Timoshenko  (6)  and  Nadai  (8)  pointed  out,  the  infinite 


series  representing  the  moments  in  the  case  of  a  concentrated  load 


do  not  converge  well.  The  numerical  computation  for  Mx,  M^,  and 


Mxy  WaS  ^one  ky  an  IBM-1620  computer.  The  computer  programs  and 


results  of  computation  are  presented  in  the  Appendix  1.  In  the 
computation,  values  of  m,  and  n  were  taken  up  to  29.  From  the 
result  of  the  calculations,  it  was  found  that  the  series  did 

3  cl 

not  converge  well  along  the  centre  lines,  x  =  -y,  y  =  — •  In  fact, 


the  series  for  M  ,  and  M  at  the  centre,  55,  were  divergent. 
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For  the  numerical  computation,  it  is  evident  from  symmetry 
that  the  computation  needs  only  be  extended  over  an  area  of  one- 
eighth  of  the  plate.  Values  of  M^,  M  ,  and  M  are  tabulated  in 
Table  2. 

Clamped  square  plate.  A  general  method  of  solution  for 
rectangular  plates  with  clamped  edges  and  an  arbitrary  loading 
has  been  developed  by  Timoshenko  (6) .  Results  of  calculation 
using  his  method  for  maximum  deflection,  moments,  and  edge  shears 
have  been  given  by  Young  (7) . 


Fig.  3.  Clamped  square  plate  with  central 
concentrated  load.  Showing  co¬ 
ordinate  axes  and  net  points. 

The  solution  of  this  particular  case  is  expressed  as 

follows : 
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in  which  w-^  is  identified  with  the  elastic  surface  of  a  rec¬ 
tangular  plate  carrying  a  concentrated  load  P  in  the  centre 
of  a  plate  and  satisfying  Navier's  boundary  conditions  w  =0 

C>2w  a2w 

and  _ W1  +  |t  1 =  0  along  the  edges.  This  first  portion  of 

d  x  by 

the  solution  contains  a  singularity  at  the  center  of  the  plate, 
where  the  bending  moments  become  infinite,  while  the  terms  w^ 
and  w^  have  no  singularities  in  the  interior  or  along  the 
edges  of  the  plate  (8).  By  using  solution  11  and  Eqs .  2,  3, 
and  4,  bending  moments,  and  twisting  moments  are  found  as 


follows : 
M 

x 


P 

2 


i  y  m  cos<px)(  ^<py)-M]^(py) j  -  y 

m=l  73 ,5  oo  m=l  ,  3 , 

cos  ( px)/^  ( py)  ( py) 


A 


M  =  TT 


P 

2 


oo 

—  )  —  cos ( fix) 

7 r  /  ,  m 


m=l ,3,5 


^  B  cos  (Py)|©1  Ox) -^l©2  Ox) 
m=l ,3,5  ^ 

d  ^  (Py)  -^(Py) 


oo 


A 


oo 


m=l ,3,5 


M 


cos  ( px) 


=  (l-tx)P 


(Py)  -1/^  (py) 


+ 


B  cos(py). 


CO 


m=l ,3,5 


trG-j^CPx)  -92  (Px) 


oo 


xy 


2 
oo 

■  I 

m=l ,3,5 


^  J  sinOx)^1Oy)  ~  y  A  sin ( px)  ^(py) 


m=l ,3,5 
B  sin(Py)  £ (Px)J 


m=l ,3,5 


where 


a  = 
m 


3  = 


m  7T 
2 

nv/r 

a 


d  = 


Poisson's  ratio 


11 


(^)(  Py)  =  (tanha. 


a 

-----  )  cosh(py)  -sinh(py)  -  tanha  sinh(By) 

cosh  a 

m 


+  m7r^  cosh(Py) 


Yi  (3y)  =  (tanham- 


a 


m 


cosh  a 


2"  )cosh(py)-2  tanha^  cosh(Py) 


m 


-tanha  sinh(py)+  sinh(py)-f  cosh(py) 

m  a  3. 

1^(Py)  =  sinh(py)-am  tanhamcosh ( Py) 

1/^(Py)=2  cosh ( py)  +  sinh(Py)-  a  tanha  cosh(Py) 

T  a.  Ill  m 

9  (f3x)=2  cosh(px)+  ■— -x-  sinh(px)--  a  tanha  cosh(px) 
1  a  '  m  m  H 

02((3x)=  ^ 


smh(Bx)-  a  tanha  cosh(Bx) 
a 

^1 (py) =( tanha^  -  — m  ■-)  )sinh(py)-  tanha^  sinh(py) 


cosh  a 


-tanha  — cosh(py)+  sinh(Py) 

m  a  v  h  j  /  a  1  ■* 

/7  „  (  P-y )  =  sinh((3y)+  cosh  (By)  -  a  tanha  sinh(py) 

<•  2  a  1  J  m  m  ^ 

j^(px)  =  sinh(px)+  cosh(px)-  a^  tanha^sinhC  px) 


A= 


ra-1 
“ 2“ 


=  B 


cosh  a 


m 


A1  =  B1  =  -0.1025  A3  =  B3  =  0.02630  A, 


Bc  =0.00420 
b 


A-.  =  B-  =  0.0015  A0  =  =  0.00055  A,  ,  =  B-,  ,=0.00021 

//  y  y  ii  ii 


A13=  B13=  °-00006 


Again  the  series  representing  the  bending  moments  did 
not  converge  well  along  the  centre  lines  x'=  0,  and  y  =  0, 

3  3  3  3 

except  at  points  (O,-^-)  /  (-j/O)  ,  (O,-^)  /  and  (-^,0)  .The  series  for 

bending  moments  were  divergent  at  the  centre  point  11.  Values 

for  M  ,  M  ,  and  M  are  given  in  Table  2.  Comuter  program  and 
x  y  xy  a 

results  of  computation  are  provided  in  Appendix  2. 

According  to  the  calculation,  the  result  seems  to  disagree 


. 


, 
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with  Timoshenko's  argument  that  the  moments  with  four  correct 
figures  may  be  obtained  by  using  seven  coefficients  A  and  seven 
coefficients  B  (6). 

Cantilever  plate.  Because  of  the  free  corners,  it  seems 
that  the  problem  of  the  bending  of  a  cantilever  plate  of  finite 
dimensions  is  difficult  to  solve  from  the  standpoint  of  a  bi¬ 
harmonic  analysis.  However,  various  approximate  methods  have 
been  applied  to  this  kind  of  problem.  The  finite  difference 
analysis  of  a  cantilever  plate  with  a  span-to~chord  ratio  of 
1:4  and  loaded  by  a  concentrated  force  at  the  middle  of  the 

longitudinal  free  edge  was  presented  by  Holl  (10).  This  problem 

_  4 

requires  a  solution  of  y  w  =  0 ,  and  all  boundary  conditions  to 
be  satisfied.  An  approximate  solution  was  attempted  by  dividing 
the  plate  into  square  nets  of  width  A  =  as  indicated  in  the 
right  half  of  the  plate  in  Fig.  4.  The  results  of  the  calcula¬ 
tions  are  given  in  Table  4.  and  represented  graphically  in  Fig. 
36,  and  Fig. 37. 


13 


Fig.  4.  Cantilever  plate  with  concentrated 
edge  load.  Showing  coordinate  axes 
and  net  points. 
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III.  EXPERIMENTAL  INVESTIGATION 

1.  The  Moire  method. 

From  the  theory  of  thin  plates,  the  equations  relating  the 

deflections  to  the  moments  are  given  by  Eqs .  2  ,  3  and  4  . 

The  moments  are  determined  by  differentiating  the  slopes.  If, 

at  some  point,  p  (x1#  y^) ,  in  a  plate  (Fig.  5)  (8),  the  moments 

M  ,  M  ,  and  M  are  desired,  it  is  possible  to  find  them  if  the 
x  y  xy 

curvatures  of  the  deflected  plate  were  known  along  the  direction 
and  yi ,  -yx . 


Fig.  5.  Determination  of  the  curvature 
and  twists  from  known  slopes. 

.  .  .  c3 

The  curvature  in  the  x  direction  is  approximately  3^*2 


.  C>w 

which  can  be  found  from  the  slope  of  the  —  versus  x  curve 

ax 
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c)  2 


Similarly  ^^2  can  ke  found  from  the  slope  of  the  versus  y 

curve.  These  curvatures,  along  with  the  constant,  D  and  jj,  can 

be  used  to  find  M  and  M  . 

x  y 

In  order  to  determine  the  twisting  moment  M  ,  it  is  neces- 

9  2.  xy 


sary  to  know  the  twist, 


w 


dx^y 


;  this  can  be  approximated  at  point 

<5w 


p,  as  shown  in  Fig.  5.,  by  plotting  the  values  of  —  along  the 
line  y-j_  -y^  and  finding  the  slope  at  point  p,  or  similarly  by 
plotting  ^  along  x^  -x^. 

Knowing  the  slopes  along  a  line  in  the  plate,  it  is  possible 

to  determine  the  deflections.  For  example,  in  Fig.  6,  if  the 

deflection  at  a,  W  is  known,  the  deflection  at  a  will  be; 

a 


aw 


w  .  =  w  + 
a  '  a 


<5w 


Fig.  6.  Determination  of  deflections. 

in  which  the  last  term  of  the  right  hand  side  is  the  area  under 

<5  w 

the  versus  x  curve  and  can  be  found  by  a  planimeter. 

Thus  once  the  slopes  have  been  found  for  all  points  in  the 

plate,  it  is  possible  to  find  deflections,  moments  and  shears 
throughout  the  plate.  The  Moire  method  provides  the  means  for 
finding  these  slopes. 

The  Moire  or  "watered  silk"  effect  is  an  optical  phenomenon 
produced  by  the  "mechanical  interference"  of  two  somewhat  simi¬ 
lar  arrays  of  dots  or  lines.  The  result  of  this  effect  is  the 

pattern  of  alternate  dark  and  white  fringes  as  shown  in  Fig.  7. 


*  ' 


. 

' 
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This  fringe  pattern  gives  the  changes  in  slopes  of  the  surface 
due  to  the  application  of  a  load,  and  from  this,  the  curvatures 
and  the  deflections  can  be  determined  as  previously  mentioned. 


Fig.  7.  Moire  effect. 

In  this  method,  the  essential  elements  are  a  ruled  screen, 
a  camera,  and  a  mirrored  model  plate. 

In  testing,  the  model  was  used  as  a  mirror  to  observe  a 
ruled  screen.  If  the  model  is  loaded,  the  image  on  the  gpund 
glass  of  the  camera  will  be  distorted.  The  change  of  the 
image  was  caused  by  the  deflection  of  the  plate  model.  The 
difference  between  the  two  images  can  be  recorded  by  superimpos 
ing  the  photographs  for  the  unloaded  and  the  loaded  state  on 
the  same  negative.  If  a  suitable  ruling  was  chosen  for  the 
screen,  "Moire"  fringes  will  appear.  With  a  special  form  of 
the  screen,  these  can  be  interpreted  as  contour  lines  for  the 
first  derivative  of  w  (1) . 

The  moments  are  liner  functions  of  the  second  derivatives 

* 

of  w.  These  are  first  derivatives  of  slope  0  and  can  easily  be 
determined  (as  illustrated  in  Fig.  5) .  For  determining  the 
moments,  the  plate  bending  rigidity  D  must  also  be  known. 

In  order  to  find  the  change  in  slope  0  of  the  model 


. 

. 
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from  the  unloaded  condition  (Fig.  8a)  to  the  loaded  condition 
(Fig.  8b) ,  the  difference  between  f  and  e  must  be  obtained, 
where  e  is  the  distance  from  the  axis  I-I  of  the  point  on  the 
surface  of  the  ruled  screen  reflected  to  the  photographic  nega¬ 
tive  by  a  point  p  on  the  unloaded  model,  and  f  is  the  distance 
to  the  point  on  the  screen  reflected  by  the  same  point  p  with 


the  model  loaded.  For  a  flat  screen, 

-r 


A=  f-< 


2L0 


1  + 


b* 


(1),  (2) 


12 


where  0is  the  change  in  slope  of  the  deflected  plate.  Expres¬ 
sion  12  will  still  be  valid  when  the  initial  slope  of  the 
model  is  not  zero  (2) . 


b)  Model  loaded. 


Fig.  8.  schematic  diagrams  for  Moire  principle. 
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If  the  point  p  on  the  plate  is  not  too  great  a  distance  from 

2 

the  axis,  the  term  ^2  will  be  very  small  and  can  be  neglected 
for  a  flat  screen.  For  example,  with  b=25  in. and  with  an  8  inch 
model  (i.e.,  maximum  c=4  in),  an  error  of  about  2.5  per  cent 
would  result  if  0were  taken  equal  to  However,  this  error 

can  be  reduced  to  a  negligible  amount  by  the  use  of  a  curved 
screen,  as  suggested  by  Ligtenberg  (1) . 

For  a  cylindrical  screen  of  radius  equal  to  3.5  Lf  the  values 
of  0  will  be  approximately  equal  to  -x— .  As  shown  in  Fig.  9,  and 
Fig.  10  (2),  if  the  model  is  of  such  a  size  that  its  extreme 
point  is  0.2  L from  the  axis,  the  maximum  error  will  be  about 
0.05  per  cent. 


Fig.  9.  Values  of  ^  at  various  points  on 
model  (for  b=L) 


* 

, 
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b  (Fig.  8)  may  not  actually  be  equal  to  l,  but, 
for  practical  reasons  may  be  greater.  The  effect  of  b  on  the 
value  of^-  was  calculated  by  Bradley  (2),  and  the  result  is 
represented  graphically  in  Fig.  10  and  tabulated  in  table  1. 

TABLE  1.  VALUES  OF  0  AND  ERROR  (E)  RESULTING 
FROM  USING  0  =  A;  FOR  CYLINDRICALLY 
CURVED  SCREEN.  RADIUS  =  3.5L. 


Y 

.1L 

.  3L 

.  5L 

.  7L 

1.0L 

0 

.49993  ^ 

.49982  A 

.49977  # 

.50035  f 

.50447  £ 

L 

b=L 

E  % 

+  .01 

+  .  04 

+  .05 

-  .07 

-  .89 

C 

.  05L 

.  151L 

.  254L 

.  363L 

.  540L 

0 

— 

.49998  ~ 

.50020  f 

i-J 

- 

.50615  f 

b=1.02L 

E  % 

— 

+  .004 

-  .04 

~ 

-  1.22 

C 

— 

.  152L 

.  257L 

.  544L 

0 

_ 

.50021  ^ 

.50053  y 

J_j 

.50858  ” 

L 

b=l ,05L 

E  % 

- 

-  .042 

-  .106 

_ 

-  .169 

C 

— 

.  155L 

.  261L 

.  551L 

0 

— 

.50056  # 

Li 

.50180  ~ 

JLi 

_ 

.51241  - 
a 

b=l . 10L 

E  % 

— 

-  0.112 

-  .359 

-  2.420 

C 

— 

.  158L 

.  266a 

— 

.  562a 

In  order  that  0  may  be  determined, A  has  to  be  measured. 
As  an  example,  in  Fig.  11a  (2) ,  the  first  exposure  was  made 
with  point  P  reflecting  a  black  line,  ,  and  if  the  resulting 
photograph  were  developed,  it  would  merely  have  recorded  a  re¬ 
flection  of  the  grid.  If  the  model  is  rotated  after  the  first 


■ 

• 

Distance  from  axis  c  Distance  from  axis 
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b=A/2L  too  small  %error  d=A/2L  too  large 

b)  Percentage  emror  resulting  from  using  d=A/2L  for 
various  values  of  b. 

Fig.  10.  Effect  of  location  of  camera  lens  on  riL/A 
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V 


T 


1  C. 


C 


a)  Unloaded  model 


b)  Loaded  model-entire 
plate  rotated  through 
angle  0 


Fig.  11.  Rotation  of  plate  to  give  photograph 
which  is  light  throughout. 


exposure  through  an  angle  0,  as  shown  in  Fig.  lib,  such  that 
a  point  p  now  reflects  a  white  line,  and  a  second  exposure 
is  made  on  the  same  negative,  the  result  will  be  a  negative 
which  has  been  sensitized  throughout,  and  a  light  print  will 


Therefore,  the  slope  0  of  the  model  can  be  found. 

If  the  model  is  rotated  through  a  greater  angle  after 
the  first  exposure  (Fig.  12)  such  that  a  dark  line  was  re¬ 
flected  to  the  negative  by  point  p,  the  film  would  be  sensi¬ 
tized  at  n  by  neither  exposure,  and  a  black  line  would  result 
at  n  on  the  final  print;  the  whole  print  would  simply  appear 
as  the  grid  itself,  with  alternate  light  and  dark  lines.  In 


,,  .  k  d 

this  case,  0  =  -  -  = 


Fig.  12.  Rotation  of  plate  to  give  photograph 
showing  original  grid. 
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Thus  whenever  the  slope  was  -j-jr,  2  (-^)  ,  ...n(— ),  the 

original  grid  would  appear  with  alternate  black  and  white  lines 


*|  ^  ^  ^ 

When  the  slope  was  (-^)  ,  (jj)  ' 
be  white. 


(—-)  ,  the  print  will 


The  model  in  Fig.  13  was  a  cantilever  plate  clamped  at 
the  top  edge.  In  the  actual  model,  the  slopes  varied  continu¬ 
ously  from  point  to  point,  and  there  resulted  on  the  prints  a 

series  of  light  regions  and  other  regions  in  which  the  original 


grid  appeared  as  shown  in  Fig.  13. 


a)  Fringe  pattern  for  slope  in  x-direction. 

b)  Fringe  pattern  for  slope  in  y-direction. 


Fig.  13.  Moire  fringes  for  cantilever  plate 

with  a  concentrated  load  at  the  mid¬ 
point  of  the  lower  edge. 

Load  =  10.019  lb.  Span-to-chord 
ratio  is  1:4.  Aluminum  sheet. 
Thickness  =  0.0630  in. 


In  these  photographs,  regions  of  the  model  for  which  0  ~ 

,  2  (-— )  ,  • . «  .  n  remained  black  which  received  the  reflection 

-  L  £  L  L 

of  a  black  line  in  both  exposures.  For  portions  which  received 

either  the  reflection  of  two  white  lines,  or  of  one  black  and 

one  white  line,  the  photograph  was  white.  For  points  on  these 
3  cj,  ,3 

lines,  0S*2  1  ~2^~2l)  ’  etc*  In  interpreting  Fig.  13  a,  since 

the  top  edge  is  fixed,  the  fringe  for  0  =  o  appears  along  the 
clamped  edge.  Moving  downward,  the  successive  dark  fringes 
show  all  points  for  which  =  0  =  ,  2(5--)  ,  etc.;  the  light 

qK  A  L  Jj 

fringes  show  points  for  which  =  0  =  (d/2L)  ,  3/2  (d/2L)  ,  and 

so  on.  For  convenience  of  plotting,  the  magnitude  of  these 
fringes  was  designated  by  order  numbers  as  shown. 

The  relative  magnitude  of  the  second  derivatives  can  be 

obtained  by  observing  the  fringe  pattern.  The  fringes  for  -~- 

2  oY 


d  w 


were  closer  together  in  tie  y  direction  where  -g— jj  was  greater. 
Similarly,  the  fringes  for  were  closer  together  in  the  x 

direction  where  was  greater. 

ax  ay 

From  the  two  photographs  of  Fig.  13,  it  was  then  possible 
to  find  the  curvatures,  and  >  an<^,  knowing  the  load,  the 

dimensions  of  the  model,  and  the  elastic  constants  of  the  materi¬ 


al,  the  moments  were  computed. 

2.  Apparatus. 

The  design  of  the  test  apparatus  was  adapted  from  reference 
(3),  For  the  purpose  of  testing  a  metal  model,  a  more  rigid 
structure  was  required.  In  this  study,  the  frame  was  built  of 
3  x  3  x  in  .steel  angles. 


They  were  welded  and  bolted  together 


, 


- 

. 
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as  shown  in  Fig.  14  and  Fig.  15.  Its  overall  size  was  72  x  36 
x  34  in.  The  distance ,L  ,  from  the  model  to  the  screen  was  25 
in.  The  distance,  b,  from  model  to  the  lens  of  the  camera  was 
26.25  in  .  i.e.,  b  =  1.05  L.  The  ruled  screen  was  36  x  36  in. 

The  grid.  The  grid  used  was  drawn  by  lining  with  India  ink 
on  a  sheet  of  white  paper.  For  obtaining  the  clearest  Moire 
fringes,  the  width  of  each  line  was  chosen  to  be  0.05  in.  and 
the  grid  spacing,  d,  was  0.10  in.  Since  the  ruling  interval 
d  must  be  the  same  everywhere,  the  ordinary  method  of  drafting 
was  believed  inadequate  to  do  such  accurate  work.  Ruling  was 
finally  achieved  with  the  help  of  a  milling  machine  and  a  par¬ 
ticular  inking  device.  A  Leroy  pen,  with  a  specially  designed 
holder,  was  held  steadily  to  a  slider  which  was  allowed  to  slide 
smoothly  on  a  steel  ruler  that  was  fixed  to  the  frame  of  the 
machine.  A  board  was  then  placed  under  the  ruler  and  was  clamped 
to  the  table  of  the  machine.  A  sheet  of  poster  paper  was  placed 
on  the  board.  After  a  line  had  been  finished,  the  table  was  then 
moved  a  distance  of  0.1  in.  ahead  by  means  of  a  micrometer  adjust¬ 
ment,  and  another  line  was  drawn,  and  so  on. 

Finally,  the  grid  was  pasted  to  the  curved  screen.  The 
radius  of  this  cylindrical  surface  was  87.5  in.  or  3.5  L.  The 
frame  of  the  screen  was  constructed  in  such  a  way  that  the  grid 
could  be  quickly  rotated  through  any  desired  angle;  index  mark- 
ings  were  placed  on  the  black  of  the  frame  so  that  the  angular 
position  could  be  easily  determined.  For  lighting  the  grid,  two 
No.  2  photoflood  bulbs  with  built-in  reflectors  were  used. 


t 

■  ■ 
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Fig.  14.  General  view  of  apparatus 


F  ig. 


15. 


General  view  of 
for  application 


apparatus.  Note  the  arrangement 
of  concentrated  load. 
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A  check  on  the  accuracy  of  the  grid  was  obtained  by  photo¬ 
graphing  the  screen  once,  and  then  after  turning  it  through  a 
small  angle  (e.g.3°)  superimposing  a  second  photograph  on  the 
first.  The  result  ought  to  be  a  pattern  of  uniformly  spaced 
straight  Moire  fringes  (1)  . 

Fig.  16  shows  that  the  screen  used  was  quite  satisfactory. 


Fig.  16.  Uniformly  spaced  straight  Moire 
fringes.  (angle  of  rotation=3°) 


The  camera.  The  camera  used  was  a  4  x  5  Crown  Graphic 
with  a  f/5.6,  180  mm  lens.  The  camera  was  rigidly  mounted 
behind  the  small  hole  (diameter=2  in)  in  the  centre  of  the 
screen.  In  order  to  avoid  any  possible  motion  during  the 
exposure  or  between  the  two  exposures,  the  shutter  was  con¬ 
trolled  by  a  metal  cable  release.  When  an  aperture  was  1/45, 
with  Kodak  Kodalith  Ortho  (type  3)  sheet  films,  good  negatives 
were  produced  with  each  of  the  two  exposures  lasting  seven 


seconds . 


. 
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Model  support.  The  plate  model  tested  was  supported  in  a 

vertical  position  in  the  test  frame.  In  this  thesis,  simply 

supported  square  plate,  clamped  square  plate,  and  cantilever 

plate  were  investigated.  For  the  clamped  square  plate,  two 
3 

•4  in.  steel  plates  were  machined  as  shown  in  Fig.  17.  The 

centre  opening  was  9*05  x  9.05  in.,  and  the  assembly  clamped 

3 

together  with  twenty  -g  in.  bolts. 

For  the  simply  supported  edge  condition,  the  fixture  was 
as  shown  in  Fig.  17,  and  Fig.  18.  The  model  was  held  verti¬ 
cally  at  the  lower  knife  edge  by  two  small  pins.  The  bolts 
were  tightened  lightly  to  draw  the  assembly  together.  On 
testing,  a  small  initial  load  was  applied  to  the  model  in  order 
to  seat  it  against  the  front  knife  edge. 

The  support  for  a  cantilever  plate  is  shown  in  Fig,  19. 

3  • 

The  plate  was  clamped  along  the  upper  edge  by  four  -g  in.  bolts. 

Loading  mechanism.  A  concentrated  load  was  applied  to  the 
model  by  means  of  a  lever  mechanism,  which  appears  in  Fig.  20 
and  Fig,  21a.  Weights  were  applied  to  one  arm  of  the  lever, 
and  the  direction  was  changed  so  that  a  horizontal  force  P  was 
applied  to  the  model.  The  two  arms  of  the  lever  were  not  equal, 
forces  were  corrected  by  a  factor  g/h. 

For  the  application  of  concentrated  load  to  the  cantilever 
plates,  a  slender  flexible  wire  was  passed  over  a  pulley,  one 
end  of  the  wire  was  connected  to  the  plate  by  a  C-clamp  as 
shown  in  Fig.  20,  and  Fig.  21b,  and  load  was  suspended  at  the 


other  end. 
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F ig.l7.Supports  for  simply  supported  square  plate 
(right)  and  clamped  square  plate(left). 
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Fig.  18.  Arrangement  for  simple  support. 


The  support  for  a  cantilever  plate  is  shown  in  Fig.  18 
It  was  secured  to  the  test  frame  by  eight  bolts. 


F  ig . 


19.  Support  for  cantilever  plate. 
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Fig.  20.  Loading  devices. 

The  model.  The  materials  used  for  this  investigation  were 
ANACONDA  cold  rolled  0,0641  in, brass  sheet  and  ALCLAD  7Q75-T6, 
QQ-A-287A  ALCOA  0.063  in.  aluminum  sheet.  In  order  to  obtain 
a  reflective  surface,  one  side  of  the  model  was  first  polished 
by  a  felt  wheel  which  was  driven  by  a  in.  hand  drill.  After 
that,  Silvo  silver  polish  was  applied.  The  plate  was  polished 
to  such  an  extent  that  the  image  of  the  grid  at  a  distance  of 
25  in.  away  could  be  seen  clearly  on  the  mirroring  surface. 

The  thickness  of  the  sheet  to  be  used  is  governed  somewhat 
by  the  fineness  of  fringe  pattern  desired  with  the  particular 
grid  being  used.  Since  small  deflection  theory  is  no  longer 
applicable  with  accuracy  if  the  deflections  exceed  about  one- 
half  of  the  plate  thickness  (2) ,  a  lower  limit  on  the  thickness 
is  imposed.  On  the  other  hand,  the  thicker  the  model,  the 


' 
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a)  Lever  system. 


b)  Pulley  system. 

Fig.  21.  Actual  view  of  loading  arrangement. 
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greater  will  be  the  loads  required  to  produce  suitable  fringe 

pattern  for  analysis.  In  this  study,  the  maximum  deflection  of 

the  plate  was  confined  to  be  less  than  ■—  in.  which  is  onehalf 

the  plate  thickness. 

3.  Calibration  of  the  material. 

It  is  possible  to  compute  the  value  of  plate  rigidity  D 

from  the  elastic  constants  of  the  model  material.  However, 

D  may  be  determined  by  using  either  a  square  model  (Fig.  22, 

and  Fig.  24)  or  a  circular  plate  (Fig.  25) .  In  Fig.  24  a  square 

plate  was  supported  at  three  corners,  and  a  concentrated  load 

was  applied  at  the  fourth  one.  The  x  and  y  axes  were  taken 

alonq  the  diagonals  as  indicated,  then  M  ,  M  ,  and  M  are 
J  J  x  y  xy 

p 

constant  everywhere:  M  =  -  M  =  -5-  >  =  0. 

x  y  xy 


Fig.  22.  Square  calibration 
plate . 


Fig.  23.  Equally  spaced  fringe 

pattern.  Load  P=0.6  lb. 
Brass  sheet  thickness 
=0.0641  in. 
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In  this  case,  the  fringes  will  be  straight  and  equidistant, 

and  spaced  a  distance  Sf  apart  (Fig.  23) 

Since;  M  =-M  =  M 

x  y 

and  from  equation  2  ,  and  3 


and  solving  for  D 
M 

D_"  (l-u)O^w 
where  the  bending  moment 


13 


and 


M  =  -  ^ 

y  2 

■\  9  'd_w  d 

dw^^y_2L 
^  y2  Sf  Sf 


Thus,  equation  (13)  becomes: 

P  S. 


D  =D=  Jf 

C  2(l-u)  {&-) 


14 


For  the  uniformly  loaded  circular  plate  with  clamped  edges, 

the  centre  deflection  has  been  given  (6)  as, 

4 

q  r 
W_  64D 

c 

where  q  is  the  intensity  of  loading,  and  r  is  the  radius  of 

the  plate.  By  measuring  the  centre  deflection,  D  for  the  cali- 

4 


bration  plate  was  found  to  be  D 


=  SUL 
:  64w. 


The  size  of 'the  square  plate  was  6  x  6  in.,  the  circular 
plate  was  9  in.  in  diameter.  Poisson’s  ratio  for  brass  was 
taken  to  be  0.33.  By  using  square  calibration  plate,  the  value 
of  Dc  for  brass  plate  was  found  to  be  468  lb-in.  (load  P=Q.61b.), 


**  Taonia 
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Fig. 24.  Actual  view  of  square  calibration  plate.  The 
plate  is  supported  at  three  corner  points  and 

a  concentrated  load  applied  at  the  fourth  con- 
er  point. 


Fig.  25.  Support  for  circular  clamped  plate.  Both 

concentrated  load  and  uniformly  distributed 
load  can  be  applied. 
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for  aluminum,  Dc=  247  (load  p=  2.377  lb.)  was  obtained  by  using 
(1=0.30.  The  value  of  Dc  (=242  lb-in.)  for  aluminum  plate  found 
by  using  a  circular  clamped  plate  with  concentrated  centre  load 
was  in  good  agreement  with  that  found  by  using  a  square  plate. 

4.  Test  procedure. 

In  order  to  ensure  that  the  deflection  of  the  plate  does 
not  exceed  one-half  the  plate  thickness,  a  deflection  -  load 
diagram  (Fig.  26)  was  first  plotted  by  loading  the  model  gradu¬ 
ally  and  measuring  the  corresponding  deflection  by  a  dial  indi¬ 
cator.  From  Fig.  26,  loads  were  chosen  so  as  to  keep  deflec¬ 
tions  of  the  model  plate  less  than  one-half  of  the  plate  thick¬ 
ness. 

In  each  test,  the  first  exposure  was  taken  with  a  small 
initial  load  on  the  model,  and  the  second  exposure  was  then 
taken  after  application  of  desired  load.  After  taking  the 
first  picture,  the  grid  was  rotated  through  90  degrees,  and 
the  fringe  pattern  for  the  other  direction  was  made.  Analysis 
was  made  upon  the  enlargements  of  these  two  negatives. 

5.  Reduction  of  data. 

In  reducing  the  data  from  the  fringe  patterns  to  determine 
the  curvatures  and  twists,  the  following  procedure  was  followed: 

Enlargements .  The  photographs  were  first  enlarged  to  full 
scale.  The  chief  quality  aimed  for  was  maximum  contrast  bet¬ 
ween  the  even  fringes  and  the  half  fringes  as  well  as  sharpness 
of  detail.  Kodak  Kodabromide,  light  weight,  A. 5,  11  x  14  in. 
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Deflection 


37 


photographic  paper  was  found  to  be  satisfactory  for  this  pur¬ 


pose. 

Preparation  of  the  picture.  Next,  grid  lines  with  spacing 
— ,  where  a  is  the  side  length  of  the  square  plate,  were  drawn 
directly  on  the  photograph,  and  then  the  fringes  were  numbered 
(Fig.  27)  .  Each  half  fringe  was  then  traced,  and  each  point 
where  it  intersected  with  a  grid  line  was  marked  with  a  very 
small  pin-hole. 


M 


versus 


Drawing  the  slope  curves.  The  curves  of  versus  x,  _ 

ax  <jx 

aw  aw 

y,-^  versus  y,  and~^  versus  x  were  next  drawn.  First  ; 


each  half  fringe  was  located  along  each  line  of  the  grid.  For 

aw 

example,  the  curves  of  versus  x  along  grid  lines  y^ ,  y^ , 

and  Y .  are  shown  in  Fig.  28,  the  curves  of  versus  y  along 

4  ax 

grid  lines  x  x, ,  x  ,  and  ,  are  shown  in  Fig.  29.  In 

o ,  l  Z  3 

drawing  these  curves,  it  was  helpful  to  draw  both  sets  simul¬ 
taneously,  for  the  results  of  one  set  was  helpful  in  locating 
points  on  the  curves  of  the  other.  For  example,  the  value  of 
at  x^  on  the  curve  y^  in  Fig.  28  is  equal  to  the  value  of 

at  y^  on  the  curve  x-^  in  Fig.  29. 


a2 


Determination  of  curvatures.  In  general,  the  curvature 


w 


was  found  by  measuring  the  slope  of  the  curves,  such  as 
in  Fig.  28,  and  the  twist  r  w  by  finding  the  slope  of  the 

axdy 

curve  as  in  Fig.  29.  These  slopes  were  approximated  at  a  point 
by  laying  a  straight  edge  such  that  it  intersected  the  curve 
at  two  points  equidistant  on  either  side  of  the  point. 
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x 


a)  For  slope  in  x-direction. 


b)  For  slope  in  y-direction. 

Fig.  27.  Fringe  patterns  for  simply  supported  square  plate 

with  concentrated  centre  load.  Initial  load=2.25  lb. 
final  load  P=14.78  lb.  Brass  sheet  thickness= . 0641  in. 


Fringe  value  r  I  unit  =  0.002 
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Fig. 28.  Plot  of  j*  ^ ...  frmge  values  against  x 

o  * 

along  y,  ,  y2  ,  y3  .  and  y4  . 


40 


Fig.  29. 


Plot  of  fringe  values 
along  x^  ,  x  {  ,  x2  ,  and  x 3  . 


against  y 
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Fig. 30.  Plot  of 


value 


against  x. 
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In  finding  the  curvature  at  the  edge,  the  values  were  ex¬ 
trapolated  by  finding  the  curvatures  at  several  points  approach¬ 
ing  the  edge  and  connecting  them  by  a  curve  (2). 

The  condition  that  bending  moments  along  the  free  edge  are 
zero  requires 

Curvature  normal  to  edge  =-(|a)  (curvature  tangent  to  the 

edge) 

Determination  of  moments.  After  the  curvatures  and  twists 
are  found,  all  data  are  available  for  the  moment  determinations. 
From  equations  2  ,  3  ,  and  4  ,  the  bending  moments  and  twist¬ 

ing  moments  are: 


M  = 

X 

-D 

(rx? + 

li 

,a2w 

.  2 
o  w. 

-D 

^B  y^  1 

u  ^  x2) 

M 

xy 

D  (1-|X) 

(d  2W  \ 
'B  aB  y 

M 

-D 

*  2 

,6  w  . 

^  2 

O  w  x 

x  = 

P 

P 

<5x2  + 

ll 

sfloi 

-D 

P 

,B  2w 

+ 

.  2 

B  wN 

If 

II 

-  D 

p 

(i-n) 

(  B  2w  ) 

p 

B  xBy 
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IV.  COMPARISONS  OF  EXPERIMENTAL  RESULTS 
WITH  ANALYTICAL  RESULTS 

1.  Test  results. 

Simply  supported  square  plate.  The  fringe  patterns  for 
a  simply  supported  square  plate  subjected  to  concentrated 
centre  load  are  shown  in  Fig.  27.  The  distribution  of  bend¬ 
ing  moments  and  twisting  moments  along  each  grid  line  is  re¬ 
presented  graphically  in  Fig.  31,  and  Fig.  32.  Numerical 
values  are  tabulated  in  Table  2. 

Clamped  square  plate.  Fig.  33  shows  the  fringe  patterns 
for  the  clamped  square  plate  with  concentrated  centre  load. 

The  distribution  of  M  ,  M  ,  and  M  is  shown  in  Fig.  34,  and 
Fig.  35.  Numerical  values  are  given  in  Table  3. 

Cantilever  plate.  Fig.  13  shows  the  fringe  patterns  for 
the  cantilever  plate  with  concentrated  load  applied  opposite 

the  clamped  edge.  Graphical  representation  of  M  ,  M  is 

x  y 

shown  in  Fig.  37.  Numerical  values  of  M^,  M^  and  w  are  given 
in  Table  4.  Also,  deflection  profiles  of  this  plate  are 
shown  in  Fig.  36. 

2.  Comparisons  of  the  results. 

Simply  supported  square  plate.  Both  experimental  and 
theoretical  results  for  the  simply  supported  square  plate  with 
concentrated  centre  load  are  given  in  the  following  table: 


« 

. 
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II  12  13  14  15  x 


Rg.3L  Distribution  of  Mx/p,  My/p  for  simply  supported  square  plate 
with  concentrated  centre  load.  P  =  l4.78lb.  Brass. 

Max.  theoretical  deflection  at  centre4  w=  — 


D 


Max.  experimental  deflection  at  cerrtre2  w= 


O.OI  1 18  P  a4 


D 
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II  12  13  14  1 5  x 
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a)  For  slope  in  x-direction. 


b)  For  slope  in  y-direction. 

Fig.  33.  Fringe  patterns  for  clamped  square  plate  with 
concentrated  centre  load.  Load  P=27.52  lb. 
Brass  sheet  thickness= . 0641  in. 
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\ 

\ 

X 

X 


Fig. 34. Distribution  of  Mx  /  P  ,  M  /  P  for  clamped  square  plate 
with  concentrated  centre  load.  P  =  27.52lb.  Brass. 


Max.  theoretical  deflection  at  centre5  w=  0.  00560P  dfyD 
Max.  experimental  deflection  at  centre5  w=0.  0052 7Pa4/ D 
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II  •  12  13  l4  15  <L 


Fig. 35.  Distribution  of  M  /P  for  clamped  square  plate  with 

xy/ 

concentrated  centre  loa d.  P  =  27.5 2  lb.  Brass. 


Deflection  in. 
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Fig.36.Defl ection  profiles  for  cantilever  plate  with  concentrated 
load  P=IO. OI9lb.  applied  at  point  31. 

Max.  theoretical  deflection  at  point  31  :  w=  0.1 8773  P  a4/  D 

Max.  experimental  de flection  at  point  31 :  w=  0.1 9040P a4/ D 
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II  12  13  14  1 5  x 


Fig.37.  Distribution  of  Mx/P,  My/P  for  cantilever  plate  with 
concentrated  load  P  =  1 0 . 0 1 9  lb .  applied  at  point  31. 
span- to—chord  r a4 i o  =  1 1  ^  ,  Aluminum. 


•  ■ 
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TABLE  2.  MOMENTS  IN  A  SIMPLY  SUPPORTED  SQUARE  PLATE  WITH 
CONCENTRATED  CENTRE  LOAD. 

POISSON  *  S  RATIO  =  0.33  (BRASS) 


NET+ 

V 

P 

M  /P 

y 

M  /P 

^ 

POINT 

Navier 

Moire 

Navier 

Moire 

Navier 

Moire 

11 

0 

0 

0 

0 

-.05840 

-.05260 

12 

0 

0 

0 

0 

-.05540 

-.04836 

13 

0 

0 

0 

0 

-.04530 

-.04497 

14 

0 

0 

0 

0 

-.02630 

-.02545 

15 

0 

0 

0 

0 

0 

0 

22 

.01140 

.01115 

.01140 

.01158 

-.05300 

-.04666 

23 

.02490 

.02783 

.02030 

.02102 

-.04460 

-.04369 

24 

.03920 

.04154 

,02470m 

.02620* 

.02837 

-.02690 

-.02757 

25 

.04610# 

.04237 

.02922 

0 

0 

33 

.04660 

.04654 

.04660 

.04527 

-.04100 

-.03818 

34 

.08000,, 

.09990* 

.07840 

.05950m 

.06240* 

.05973 

-.02900 

-.02969 

35 

.10626 

.07232 

0 

0 

44 

.11800 

.11054 

.11800n 

.13000* 

.11266 

-.03000 

-.03182  ! 

45 

.17700 

.18292 

.14388 

0 

0 

55 

* 

.31922 

* 

.31922 

0 

0 

Observing  this  table  the  maximum  error  is  approximately  16 
per  cent  at  net  point  35  for  M  .  The  maximum  deflection  at 

2  y 

Pa  ^ 

centre  is  wmax  =  0.01118  ■  ■  ,  approximately  4  per  cent  less  than 

theoretical  result. 

Clamped  square  plate.  Experimental  results  were  compared 
with  Young's  result  (7).  The  result  of  calculation  is  given 
in  Table  3. 


* 

# 

+ 


Series  diverges  at  this  point. 

Series  converges  slowly  at  this  point 
For  net  points  see  Fig.  2. 


. 

■ 

' 
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TABLE  3.  MOMENTS  IN  A  CLAMPED  SQUARE  PLATE  WITH  CONCENTRATED 
CENTRE  LOAD.  POISSON'S  RATIO  =  0.33 
MATERIAL:  BRASS, 


NET+ 

Mx/P 

M  /P 

y/ 

M  /P 
xy' 

POINTS 

Young 

Moire 

Young 

Moire 

Young 

Moir& 

11 

* 

.31665 

* 

.31665 

0 

0 

12 

# 

.07364 

# 

.12125 

0 

0 

13 

# 

-.00204 

# 

.04420 

0 

0 

14 

-.05100 

-.05077 

-.00200 

-.00221 

0 

0 

15 

-.12500 

-.10986 

-.04100 

-.04622 

0 

0 

22 

.06300 

.05816 

.06300 

.05986 

-.02910 

-.02506 

23 

# 

.00374 

.02976 

.02976 

-.02460 

-.02210 

24 

-.04500 

-.04677 

-.00540 

-.00476 

-.01630 

-.01481 

25 

-.10500 

-.09098 

-.03400 

-.02993 

0 

0 

33 

-.00196 

-.00181 

-.00196 

-.00181 

-.03016 

-.02688 

34 

-.03000 

-.03163 

-.01000 

-.01037 

-.02237 

-.02005 

35 

-.05700 

-.05782 

-.01890 

-.01900 

0 

0 

44 

-.01480 

-.01497 

-.01480 

-.01497 

-.01637 

-.01731 

45 

-.01080 

-.11900 

-.00360 

-.00393 

0 

0 

55 

0 

0 

0 

0 

0 

0 

By  comparing  the  results  from  table  3,  the  maximum  error 
is  also  16  per  cent  at  net  point  25  for  M  .  Maximum  deflection 

y  2 

pa^ 

at  the  centre  of  the  plate  is  w  =  0.00527  ■  •  ,  it  yields  6 

x  max  D  u 

per  cent  less  that  obtained  by  Young. 

Cantilever  plate.  The  finite  cantilever  plate  subjected 
to  a  transverse  loading  is  one  problem  for  which  an  exact  solu¬ 
tion  has  not  been  achieved.  However,  Holl's  (10)  finite-differ¬ 
ence  solution  provides  a  check  for  this  case. 


*  Series  diverges  at  this  point. 

#  Series  converges  slowly  at  this  point 
+  For  net  points  see  Fig.  3. 
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Regarding  the  accuracy  obtainable,  it  is  observed  that 
maximum  deviation  is  about  20  per  cent  at  net  point  2  l  for 
M  ,  and  18  per  cent  at  net  point  11  for  w.  Values  of  bend- 
ing  moments  and  deflections  are  given  in  Table  4. 

In  comparing  the  results,  it  is  to  be  noted  that  Holl's 


solution  is  an  approximate  solution. 

TABLE  4.  MOMENTS  AND  DEFLECTIONS  FOR  A  CANTILEVER  PLATE 
WITH  A  CONCENTRATED  EDGE  LOAD  AT  THE  MID-POINT 
OF  THE  FREE  EDGE.  POISSON'S  RATIO  =0.30 
MATERIAL:  ALUMINUM.  SPAN-CHORD-RATIO  =1:4 


NET+ 

Vp 

M  /P 

»/£§? 

POINTS 

Ho  11 

Moire 

Ho  11 

Moire 

Ho  11 

Moire 

11 

-.14902 

-.14640 

-.49672 

-.48807 

0 

0 

12 

-.11206 

-.11512 

-.37352 

-.38405 

0 

0 

13 

-.06802 

-.07149 

-.22672 

-.23861 

0 

0 

14 

-.03780 

-.04067 

-.12600 

-.13558 

0 

0 

15 

-.01498 

-.01430 

-.04992 

-.04733 

0 

0 

21 

+.04696 

+.05620 

-.21714 

-.21643 

.01006 

.00996 

22 

-.03931 

-.03352 

-.16670 

-.17156 

.00756 

.00843 

23 

-.05539 

-.05423 

-.11477 

-.10600 

.00459 

.00478 

24 

-.03444 

-.03796 

-.07780 

-.07592 

.00255 

.00276 

25 

0 

0 

-.06432 

-.06089 

. 00101 

.00119 

31 

+.37703 

+.39711 

0 

0 

.03041 

.03090 

32 

+.00186 

+.00197 

0 

0 

.022022 

.02360 

33 

-.06766 

-.07420 

0 

0 

.01355 

.01473 

34 

-.05070 

-.04856 

0 

0 

.00809 

.00860 

35 

0 

0 

0 

0 

.00488 

.00548 

Pa' 

D 


The  maximum  experimental  deflection  is  w  =  0.1904 

^  max 

which  is  1.59  per  cent  greater  than  that  given  by  Holl's 


solution. 


+  For  net  points  see  Fig.  4 
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V.  DISCUSSIONS  AND  CONCLUSIONS 

The  use  of  plastic  materials  for  the  Moire  Method  was 
well  described  by  Ligtenberg  and  Bradley.  One  of  the  dis¬ 
advantages  of  using  plastic  materials  is  that  the  modulus 
of  elasticity  varies  with  humidity  and  temperature.  This 
will  be  an  obstacle  in  applying  Moire  method  to  some  other 
plate  problems,  for  example,  thermal  stresses  in  plates. 

Since  the  values  of  moduli  of  elasticity  for  metals  may  be 
considered  as  constant  for  practical  purposes,  it  is  pos¬ 
sible  to  take  advantage  of  this  property.  And  once  the 
material  is  calibrated,  there  is  no  need  to  repeat  calibrat¬ 
ing  for  the  same  kind  of  material  in  other  test.  In  such 
case,  much  labour  can  be  saved  during  each  test.  Therefore, 
the  use  of  metal  models  is  also  particularly  attractive 
from  the  standpoint  of  saving  time. 

i 

In  applying  the  Moire  method,  several  possible  sources 
of  errors  may  be  listed  as  follows: 

a.  Membrane  action  of  the  plate  due  to  large  deflection. 
According  to  Sturm  and  Moore's  experimental  result 
(9) ,  for  plates  in  which  the  deflection  does  not 
exceed  one  half  the  thickness  of  the  plate,  less 
than  5  per  cent  of  the  load  may  be  assumed  to  be 
carried  by  membrane  force. 

b.  Due  to  calibration.  The  value  of  D  found  by  using  a 

square  calibration  plate  with  a  concentrated  load 

p  =  2.3771b -at  one  corner  was  D  =  247  lb-in.;  by 

c 


‘ 

' 

mo 
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using  a  circular  clamped  plate  with  a  concentrated  load 

P=  14  lb.  at  the  centre  of  the  plate,  Dc  =  242  lb-in? 

for  a  circular  clamped  plate  with  uniformly  distributed 

lb 

load  q  =  1.083  -±-~ D  =  239  lb-in.  The  maximum  error 

•  2c 
in 

is  approximately  5  per  cent. 

c.  Due  to  the  use  of  the  incorrect  value  of  Poisson's  ratio 
for  the  particular  plate  tested.  For  example,  if  |i  is 
taken  as  0.32  instead  of  0.33,  the  maximum  difference 
for  moments  by  using  \±  —  0.3  2  was  11  per  cent  greater 
than  using  \x  =  0.33. 

d.  Due  to  experimental  errors  in  the  measurement  of  load-- 
it  is  estimated  as  about  1  per  cent. 

e.  Due  to  reduction  of  data--it  may  be  as  great  as  4  per 
cent . 

f.  Due  to  possible  imperfect  realization  of  a  boundary 
condition. 

g.  Due  to  shear  deflection  not  accounted  for  in  the  theory. 

As  Ligtenberg  (1)  pointed  cut,  accurate  results  are  impos¬ 
sible  near  a  point  load.  This  is  caused  in  part  by  the  fact 
that  the  theory  is  applicable  to  the  middle  plane  of  the  plate, 
whereas  measurements  are  taken  from  the  surface.  Furthermore, 
even  the  elementary  theory  is  not  applicable  near  a  concentrat¬ 
ed  load  because  the  stresses  normal  to  the  plane  of  the  plate 
cannot  be  assumed  to  be  zero.  Observing  table  2,  and  3,  the 
series  for  bending  moments  do  not  converge  well  along  the 
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centre  line  of  the  simply  supported  square  plate  and  the 
clamped  square  plate,  and  diverge  at  the  centre  of  each  plate. 

Shearing  forces  and  are  obtained  from  the  third 

derivatives  of  the  deflections.  It  was  found  that  the  values 
.3  N  3 

n  W  n  W 

of  ,  and  cannot  be  measured  with  sufficient  accuracy 

by  the  Moire  method. 

In  comparing  experimental  results  with  theoretical  results, 
it  is  observed  that  the  maximum  error  is  16  per  cent  for  the 
bending  moment,  and  14  per  cent  for  the  twisting  moment  in  the 
cases  of  simply  supported  square  plate  and  the  clamped  square 
plate;  for  the  cantilever  plate,  20  per  cent  deviation  was  ob¬ 
served  for  M^,  about  10  per  cent  for  M^,  and  about  18  per  cent 
for  deflection.  It  is  of  interest  to  note  that  the  finite 
elements  are  used  in  the  formation  of  difference  equation  net¬ 
works.  Consequently,  the  results  obtained  are  only  approxim¬ 
ate.  This  approximation  can  be  improved  as  the  elements  become 
smaller.  For  the  case  of  cantilever  plate,  the  accuracy  cannot 
be  stated  exactly. 

Moire  method  yields  approximate  solution  to  laterally  loaded 
plate  problems  with  complicated  boundary  conditions.  From  the 
viewpoint  of  engineering  application,  the  Moire  method  is  be¬ 
lieved  to  be  a  useful  tool  for  obtaining  valuable  informations 
concerning  bending  moments,  twisting  moments,  and  deflections 
in  bending  of  plates  with  complicated  boundary  conditions. 

Among  the  problems  of  bending  of  plates,  those  which  might 

I 

be  studied  effectively  by  the  Moire  method  are: 


. 
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a.  Thermal  stresses  in  plates--since  the  modulus  of  elas¬ 
ticity  may  be  assumed  as  constant  within  certain  range 
of  temperature,  it  is  possible  to  take  advantage  of 
this  particular  property  to  apply  Moire  method  to  therm¬ 
al  stress  problems  in  thin  plates.  Preliminary  inves- 

#  ,  i 

tigation  shows  that  Moire  method  can  be  applied  to 
thermal  stress  analysis. 

b.  Moments  and  deflections  in  some  elastically  supported 
plate  problems. 

c.  Study  of  plates  with  stiffeners. 

d.  Plates  with  variable  thickness. 

i 

It  is  believed  that  the  Moire  method  would  offer  useful  re¬ 


sults  in  these  studies  and  would  provide  valuable  informations 
for  application  in  practical  designs. 
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APPENDIX.  COMPUTER  PROGRAMS  AND  DATA 


1.  Simply  supported  square  plate  with 
concentrated  centre  load. 


•  • 


I  MECH.  ENG.  H.  Y.  TSO  923114 

DIMENSION  A(20,20)*B(20»20)*C(20,20> 

1  READ  2  *  I  9  IX, JY,U 

2  FORMAT ( 3  I  5 *F1 0. 5 ) 

PUNCH  5 

3  FORMA  T ( 2X  2HMN  *  4X 1 HL , 4X 1 HK , 1 0 X3HT XM , 1 7X3HT YM , 1 7 X4HT X YM / / ) 
DO  30  K  =  1  »  J  Y 

DO  30  L=K  *  IX 
DO  10  TM  = 1 ,  I 
M  = 2*  \  m- i 
AM  =M 
BL=L-1 

PSM=S INF ( AM*1 . 570795 ) 

RUM=AM*BL*0. 392ARQ 
QSM  =  S  INF ( RUM ) 
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no  10  I N= 1  »  I 

A  N  =  N 
CK=K-1 
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RUN  =  AN*CK*0  . 392 A 90 
QSN=S INF ( RUN ) 
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2.  Clamped  square  plate  with  concentrated  centre  load. 


..I  f-T  H  •  F  NG  •  H  •  r  .  ISO  923] 14A 

•  •LOAD  FOPGO  r LOCK  1  ?u0 

Df  MENS  ION  AX  < 10 )  *8 X ( 1  )  t CX ( 10 ) > AY (I  ) » B  ( 1  )  ,CY ( 10 ) 
DIMENSION  Axn 1 v) >HXr ( lu) »Cai (lG)»A(lu)*Bilo) 

1  READ 

2  FORMAT ( 31 5* FI  0,5) 

PUNCH  9 

5  FORMA  i  (  2  a  2  H  I  M4  a  1 H L  »  4  A  1 H i\  »  1  0  A  3  H  xM,i7x3H  .  M*17;-4H.;,  I//) 

DO  4  J= 1 «  I 

4  READ  ? 9 A (  J )  , R  (  J  ) 

3  FORMA  I  ( 2F  1  Q • 6 ) 

DO  20  K=1,JY 

DC  20  L=K,JX 

DO  1 0  mm= 1 , I 

M= P*MM- 1 

C  M  =  M 
DL=L~1 
D  K  =  K  -  1 

PM L=CM*DL *0.392699 
PMK=CM*DK*0 . 392699 
COSL=COSE ( PML ) 

COSK=COSF ( PMK ) 

S I NL  =  S T  NF ( PML  ) 

S I NK=STNF ( P^K ) 

A  P  =  CM*  1  .E7J-7QK 
AN  =  — AP 
EXP=EXPF ( Ah ) 

EXN=EXPF( AN) 

PuP  =  E aP  +  E  aN 
PuN=FaP-EaN 
COSHM=PJP/2 • 0 
SI N  H  M  =  P  u  N / 2 • U 
i A N H M  =  INHM/CQ  mM 
PNL =-PM( 

PNK=-PMK 
FXPK=RXPF ( PMK ) 

FXNK=FXPF ( PNK ) 

EXPL=FXPF ( PML ) 

E  aNL  =  E a PR  i PNL ) 

YuP  =  t  aPK+E  aNix 
YUN=EXPK-EXNK 
XUP=EXPL+EXNL 
XUN=EXPL-EXNL 
COSHY=YUP/2.0 
S I NHY  =  YUN/ 2*0 
CO  :,Ha=  A  p  /  2  . 

S I NHa=XuN/2. ^ 

H  j  = AP* i ANHM 
i CA  =  j ANHM-AP/ CO^HM**? 

OSHy  =  PMK*S I NH Y 
TH y =QSHy-Hl *C05HY 
OoHa=PML*  INHa 
i  Ha  =  QoHa~H  ,  *CO..H,v 
N= ( M-l ) /2 
S I = ( - 1 ) **N 
RCOSH-1 • u /COSHM 
A  M  =  5  I  *  A (MM ) *RC  OoH 
RM  =  S I*R ( MM ) *RrQSH 
VSHY= i ANHM*QSH - 
OCHY=PMK*COSH Y 
i A  r  = , CA*COoH t 
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(  /  .1 _ .  (  t  _ _ Ll  t i L_  Ah . -1  .  .HM. Ill _ 
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(•i*l  )  '  £  Tyi  Q 

_ _  ?  H^flUq _ _ 

I,  •  -  .  Hr  [  (  /  H  f  *  JH  I  A  !  H\  A  )  '  AMSU  i 

   . .  I  >  r  -L  a  •  * 
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_  (  ■>.  r  K  )  r  AMOQH  r _ 

vt*  f-y 

•'I  »  '  ~  1  V  '■» _ . 

T  ,  r  sWM 

r-vM*<;=M 
V  =  M  ^ 

_ r-  J=  ,)'  < _ 

r  -  n  =  >  f\ 

_ V-'  .  Jr  '  :  •  •■sJi'H _ 

■  C  ^  P  „  *  l#t  sNf-’M 

_ '  )  -  I  _ 

(  XMH  )  1?; 

_ _ (  J^h  )  ^  f  ''  =  JK  ]  J. _ _ 

(  >^s  )  i  ? -=y\/\  i  ? 
aoryr?,  r»M->~qA 

(  />  A.  IiMa  =  H  A  "" 

^lA  1  +  MX  l=MuM 

,  ,q\M  q  =  MH'n-> 

jMq  — =  j^q 

>Mq_-  Nj/q _ 

( >wq ) 

_ ■ _ f  !  "  ^  )  1  1 -I  <  "  =  ■-"  ‘  ^ _ 

(  JMM  )  iqxis  iq/q 

(  J Y' h  )  jqx  -1  =  J/X3 

JflX  -l+JSX  1  =  HUX 

_ lyix  H-  ) g X ^  M U X _ 

.S\q>  Y  =  YH .  0~) 

_ y  =  YKi^  1  ? _ ; 

,s\q  =ah  cn 

_ v!'  K  =  AH/I _ 

MHtfA,*qA*  H 

■>  .  qA-V-M/*  A  ,  sA  ; _ . 

V  H !/  I  *  'JMq-yH^O 

_ ,H  rVH|H-fH?P=yHT _ 

xHM  *j!  q*xH  o 
< r-w) *n 


Hi  =  i  A  i  —  o  I  NH  .  .  H  +  0  C  H 

FHXX-2 • 0*CO  HX+Q  HX-H1*C0  Ha 
CH  f=  i  ANHM*CO~H i 

i  H  r  Y2-  '  AY  — 2«  Q*C HT  — V  >H  ,  +..  I  NH  ,  +Q(  H  * 
I"  H  r  =  2  •  U  "  Cu.H  i+O  H  i  -  A  k  «-  C  H  i 
H  i  =  i  ANHM*o  1  NH  . 

[' H  X  Y  1  =  l  CA*o  I  NH  i  —  bn  ■  “Cm  i  *P  Mi\  +  u  ; ; 

Ha;  2=  I  NH  /  +QCH  -  H  *AP 
Of HX=PML*C0SHX 
FHX=HT  *S  I  NHX 
1  h  X  r  A  =  S  T  N H a  +  Q r H A  —  Ha 
A X ( MM ) =  CO  L *( H  <  —  u * ; H ,  2 ) /  C M 

RX ( MM ) =AM*COoL* l i H  -u* i H i i ) 

CX ( MM ) =RM*COGK* l I HAa-  *  H.J 
AY ( MM ) =COoL* ( y - , H  , 2 ) /CM 
By (MM) =AM*COoL*{ o* . H , -  H  ) 

C i ( MM ) =RM*C0SK* ( J* i Ha a- . Ha ) 

Aa  i  (  MM  )  =o  I  NL*  i  Ha  ;  1  /  CM 
B a  Y ( MM ) =  AM*  I NL*  H X t x 
C X  Y  (  M M  )  =  B M ,  I  N  L  *  .  H  A  3 
S  U  X  M  =  • 

SU Yw  =  0 . 0 
S  U  X  Y  M  =  0  •  0 
I  M  =  1 
R  A  A  =  ./  • 
b  B  X  =  u  •  u 

oC  A: j  •  v 

b  A  i  =  •  u 

S  B  y  =  u  •  ^ 

S  C  Y  = 

3 A  a i =o*0 
S  B  X  Y  =  0  •  0 
S  C  X  =  <  •  0 
00  V;  MM=  1  *  I  M 
SAX=SAX+AX(MM) 

SBX=SBX+BX ( MM ) 
oCX=oCX+CX( MM ) 

S A  <  =  jA i +A  y ( MM ) 

SB  i  =  oB r  +BY  ( MM ) 

SO  =SC  t  +C  i  l MM ) 
r  A  a  i  =o A a , +Aa  r  ( MM ) 

S  B  X  Y  =  S  B X Y  +BXY ( MM ) 

S  C  X  Y  =  SC  X Y  +CX  Y ( MM ) 

P= 1.0/3. 141*9 
„  L  X M  =  P *  -  A  X  — ..  B X  +  CX 
SUYM=P*SAY-SBY+ dCY 
SUXYM=P*SAXY-SP XY-SC XY 
G=  0 . 5 

GU  = (  1 . 0— U ) *0 . 3 
T  XM  =  G*SUXM 
T Y  M  =  G*SU  Y  M 
T a i M=Gu*3UX ! M 

Punch  6* im»l>k> i xm*  m > ,  m 

FORMAT (315  *3t2G.8 ) 

IF (  I  M  —  I  )  14*20*20 
IM  =  IM+1 
T  X  M  =  0 . 0 
T  Y  M  =  0  •  0 
TX YM  =  n#  0 
SU  XM  =  0 . o 
SUYM=G . 0 
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_ <  n;m  >  x3+OcgX7>a _ 
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O  A  I  M  —  •  v 

GO  TO  15 
2  ^  CONTINUE 
GO  TO  1 
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